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Abstract 

We prove Lipschitz stability in the determination of the scalar potential in 
the dynamic Schrodinger equation in an infinite waveguide from one boundary 
observation and one (arbitrarily small with respect to the infinite direction of the 
waveguide) internal measurement of the solution. 
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1 Introduction 

1.1 What we are aiming for 

In the present paper we consider an infinite waveguide = uxl, where u is a connected 
bounded open subset of R™" 1 , n > 2, with C 2 -boundary du. Given T > we consider 
the following initial boundary value problem 

-iu' - Au + q(x)u = 0, in Q = (0, T) x Q, 

u(0,x) = Uo(x), iGfl, (1.1) 
u(t,x) = g(t,x), (t,x) G £ = (0,T) x r, 

where r := du x R and the sign ' stands for J^. Here «o (resp. g) is the initial (resp. 
boundary) condition associated to (1.1) and q is a function of x G Q only. 

Since T is unbounded it is worth making the boundary condition in the last line of 
(1.1) more precise. Writing x = (x', x n ) with x' = (xi, . . . , x n .-i) G oj for every x G Q 
we extend the mapping 

C °°((0, T) x R; H 2 (a;)) — )• L 2 ((0, T) x R; H 3/2 (du))) 

v ^ [(t, x n ) e(0,T)xB4 v{t, -, x n )\du], (1.2) 
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to a bounded operator from L 2 ((0, T) x E; H 2 (w)) into L 2 ((0, T) x E; H 3 / 2 (<9w)), denoted 
by 7o. Then for every w G C°([0, T}; H 2 (f2)) the above mentioned boundary condition 
reads 70 w = g. 

The main purpose of this paper is to prove stability in the determination of the 
scalar potential q from both boundary and arbitrarily small (with respect to the infinite 
direction x n of the waveguide) internal measurements of the solution u to (1.1). 

1.2 Published papers 

The problem of stability in determining the time-independent electric potential in a 
Schrodinger equation from a single boundary measurement was treated by Baudouin and 
Puel in [1]. This result was improved by Mercado, Osses and Rosier in [8]. In these two 
papers, the main assumption is that the part of the boundary where the measurement 
is made satisfies a geometric condition related to geometric optics condition insuring 
observability. This geometric condition was relaxed in [2] under the assumption that the 
potential is known near the boundary. In all the above mentioned papers the Schrodinger 
equation is defined in a bounded domain. 

In this paper we investigate the problem of determining the scalar potential in the 
Schrodinger equation defined in a closed waveguide, which is an unbounded domain. It 
turns out that there is a very small number of articles dealing with inverse boundary 
value problems in an unbounded domain in the mathematical literature. In [7] Li and 
Uhlmann prove uniqueness in the determination of the scalar potential in an infinite 
slab from partial Dirichlet-to-Neumann map. In [4] Cardoulis, Cristofol and Gaitan 
obtain Lipschitz stability from a single lateral measurement performed on one side of an 
unbounded strip. For an inverse boundary value problem stated in a waveguide geometry 
we refer to [5, 6] where stability is claimed for various coefficients of the Schrodinger 
equation from the associated Dirichlet-to-Neumann map. Here we investigate the same 
type of problems in absence of any information given by the Dirichlet-to-Neumann map. 

1.3 Main results and contents 

In this section we state the main result of this article and briefly comment on it. 

Theorem 1.1. For M > 0, i > and a > fixed, let u E H 4 (fi; E) obey 

u (x) >a>0, x e Q, (1.3) 

let qj E Q M := {q E ^ 2 '°°(ft;E), \\q\\ w ^(n) < M}, j = 1,2, fulfill 

qi{x) = q 2 {x), x e u x (E\ (-£,£)), (1.4) 

and letuj denote the C\[0, T]; H 2 (fi) n Hj(ft)) D C 2 ([0, T}; L 2 (tt)) -solution to (1.1) as- 
sociated to uq, g = 70G and qj, where 

G(t,x) := ito(x) +it(A - q 2 )u (x), (t,x)eQ. (1.5) 
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Then for every L > i, there exist a subboundary T* C du x (-L, L) and a constant 
C > depending only on L,T,M,u and T*, such that we have 



Ql ~ <?2||l 2 (C) < C (\\dvill! - «2)||H 1 (0,T;L 2 (r,)) + |K ~ U 2 || H l (0,T;Hi (wx &£))) ) ■ C 1 - 6 ) 



Under the prescribed conditions (1.3)— (1.5), Theorem 1.1 claims Lipschitz stability in 
the determination of the scalar potential appearing in the dynamic Schrodinger equation 
in Q from two different observations of the solution u to (1.1). The first one is a lateral 
measurement of d v u on some subboundary of duo x (— L, L). The second observation is 
an internal measurement of u which is performed in each of the two "slices" S-(L) : = 
cu x (—L,—£) and S + (L) := oj x (£,L) of Q. Although the M n -Lebesgue measure of 
S±(L) can be made arbitrarily small by taking L sufficiently close to £, this observation 
cannot be removed from the rhs of (1.6) since the prefactor (ie the constant C) tends to 
infinity as L goes to I. The occurence of this internal observation of u in (1.6) arises from 
the unbounded geometry of Q. More precisely, the derivation of the stability inequality 
(1.6) being by means of a global Carleman estimate for the Schrodinger equation in a 
bounded domain, the strategy used in this paper involves a cut off function with first 
derivative supported in (-L, —£) U (£,L), which gives rise to the measurement of u in 



The paper is organized as follows. Section 2 deals with the direct problem associ- 
ated to (1.1). Namely §2.1 gathers existence and uniqueness results for the solution 
to the dynamic Schrodinger equation in the infinite domain Q and §2.2 is devoted to 
the study of the direct problem for the linearized system associated to (1.1). In §2.3 
the corresponding solution is suitably extended to a function of [— T, T] x Q which is 
continuous wrt the time variable t. This is required by the method used in the proof 
of the stability inequality (1.6), which is the purpose of Section 3. It is by means of a 
Carleman estimate for the Schrodinger equation stated in §3.1. Finally §3.2 contains 
the completion of the proof of Theorem 1.1. 

2 Analysis of the direct problem 
2.1 Existence and uniqueness results 

This subsection gathers two existence and uniqueness results needed for the analysis of 
the direct problem. We start by recalling from [6] [Proposition 2.1] the following: 

Proposition 2.1. Let M > 0. Then for all q E Q M , v £ Hj(ft) H H 2 (fi) and 
f £ W^^T-L 2 ^)), there is a unique solution v £ Z = C([0, T}; Hj(fi) fl H 2 (f2)) n 
C 1 ([0, T]; L 2 (f2)) to the boundary value problem 



S ± (L). 
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Moreover we have 

IMUo < C (\\v \\ H 2 {n) + ||/||wi>l(o,T;L 2 (n))) , 

for some constant C > depending only on uj, T and M . 
Set 

^o:=7o(^ 2 ' 2 (0,T;H 2 (fi))). (2.7) 
The space J^o is equipped with the following quotient norm 

IMUo = mf{||G|| w2 , 2(0 ,T;H^)); G G l^ 2 ' 2 (0,T; H 2 (fi)) satisfies lQ G = g}, 
in such a way that every g G %o admits an extension G G W^ 2 ' 2 (0,T; H 2 (f2)) obeying 

||G1|w 2 > 2 (0,T;H 2 (n)) < 2 ||^||s" . 

Further, put 

& ■= {K g) G H 2 (fi) x Jfo; u = g{0, •) on T}. (2.8) 

Then, setting u — v + G, where v is defined by Proposition 2.1 for / = idtG + AG — VG 
and v o = uq — G(0, .), we obtain the following existence and uniqueness result, which is 
similar to [6] [Corollary 2.1]. 

Theorem 2.2. Let M > 0. Then for every q G Qm and (uo,g) G ££ there is a unique 
solution u G Z = C([0, T]; H 2 (f2)) fl C^QO, T]; L 2 (f2)) to £/ie boundary value problem 
(1.1). Moreover the estimate 

\\u\\ z < C(||wo|| H 2(n) + IMU Q ), ( 2 - 9 ) 
holds for some some positive constant C depending only on u, T and M. 

2.2 The linearized problem 

Let M, qj, j = 1, 2, u , # and G bethe same as in Theorem 1.1. Since G G ty 2 ' 2 (0, T; H 2 (Q)) 
and g(0, •) = u on T, we have (uq, g) G Jzf . Applying Theorem 2.2 there is thus a unique 
solution Uj G Z,j = 1,2, to the system 

{— iu'j — Auj + qj(x)uj = 0, in Q, 

Uj(0,x) = u (x), ieO, (2.10) 
Uj(t,x) = g(t,x), (t,x) G E. 

Further, differentiating (2.10) with respect to t for j = 2, we obtain that u' 2 is solution 
to the following boundary value problem 

—iu 2 ' — Au' 2 + q2(x)u 2 = 0, in Q, 

1*2(0, x) = u (x), xeQ, (2.11) 
u' 2 (t,x) = g'(t,x), (t,x)eS, 
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where z/ = i(A — q 2 )uQ. Since G'(t,x) = Vo(x) for a.e. (t, x) G Q by (1.5) then it is 
true that G" G W 2,2 (0, T; H 2 (f2)) and consequently (vo,g') G Therefore m' 2 G Z from 
(2.11) and Theorem 2.2. As a consequence we have 

u 2 G C^MjH 2 ^)) nC 2 ([0,T];L 2 (fi)). 

Moreover, it follows from (2.10) that u = u\ — u 2 is solution to the linearized system 

—iv! — Au + q\u = f, in Q, 

u(0,x) = 0, (2.12) 
u(t, x) = 0, (t, x) G E, 

with / = (q 2 - qi )u 2 . Since u 2 E Z and gi - q 2 G L°°(fi) then / G W^'^T; L 2 (fi)) so 
that we have u G Zo by Proposition 2.1. 

Last, we deduce from (2.12) that v = u' satisfies 

-iv' — Av + q±v = /', in Q, 

v(0,a;) = i(g 2 — qi)(x)u (x), xEVt. (2-13) 
f (t, x) = 0, (t, x) G S, 

with /' = (q 2 — q\)u' 2 G Vr 1 ' 1 (0,T; L 2 (f2)). From the identity q\ = q 2 on T we get that 
i(q 2 — qi)uo G Hj(f2) n H 2 (f2). Therefore t> = — w 2 G Z by Proposition 2.1. Bearing 
in mind that u' 2 G Z we deduce from this that m' x G Z. Summing up, we have obtained 
that 

Uj G C l ([0, T}; H 2 (fi)) n C 2 ([0, T]; L 2 (fi)), j = 1, 2. 



2.3 Time symmetrization 

In this subsection we prove that v may be extended to a solution of (2.13) over the time 
span [— T, T] . 

We extend v (resp. /') on (— T, 0) x Q by setting v(t,x) = —v(—t,x) (resp. f'(t,x) = 
—f'(—t,x)) for (t,x) G (— T, 0) x f2. Since the initial conditions w(0, •) and /'(0, •) are 
purely complex valued according to (2.13), the mappings t \- > v(t,x) and t n- f'(t,x) 
are thus continuous at t = for a.e. 16SI. Therefore, we have w G C°([— T, T]; Hj(f2) fl 
H 2 (fi)) fl C^Q— T, T]; L 2 (f2)) and /' G W^-T, T; L 2 (fi)). Moreover v is solution to 
(2.13) over the whole time span (— T, T): 

-iv' -Av + q x v = f, in (-T, T) x Q, 

v(0,x) = i(q 2 — qi)(x)u (x), iGSl. (2-14) 
v(t,x) = 0, (t,x) G (— T, T) x 
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3 Stability inequality 
3.1 Global Carleman estimate 

In this subsection we recall a global Carleman estimate for the Schrodinger operator, 
which is borrowed from [1] [proposition 3]. 

Given an arbitrary bounded domain O C W l with C 2 boundary, we consider the 
Schrodinger operator L acting in (C^°) / ((— T, T) x O), 

L := -id t - A + q, (3.15) 

with scalar potential q G L°°(0). 

Next we introduce an open subset T of dO, and a function (3 G C 4 (0; R+) satisfying 
the following conditions: 

Assumption 3.1. 

(%) 3C*o > such that the estimate |V/3(x)| > Co holdsfor all x G O; 

(%%) d u ]3(x) := Vj3{x).v{x) < for all x G dO\T ; 

(m) 3Ai > 0, 3e > such that we have A| Vf3(x).(\ 2 + D 2 (3((, () > e\(\ 2 for all ( G W 1 
and A > Ai. 

There are actual functions (3 fulfilling Assumption 3.1, such as O 9 i i- > \x — Xq\ 2 
where xq G M. n \ Q is fixed. 
Further we put 

(3 := (3 + K, where K := Tnll/JHoo for some m > 1, (3.16) 
and define the two following weight functions for A > 0: 

e \/3(x) 2\K _ A/30) 

(p(t,x) = — — r and rj(t,x) = — — -, (t,x) G (— T, T) x O. (3.17) 

ry 1 (T + t)(T-t) /v ; (T + t){T-t) y 1 y ' K ' 

Finally, for all s > we introduce two operators acting in (C X) )'((— T, T) x O), 

Mi := id t + A + s 2 \ Vr]\ 2 and M 2 := isrj' + 2sVr/.V + s(Arj). (3.18) 

Mx (resp. M 2 ) is the adjoint (resp. skew-adjoint) part of the operator e~ sr, Le sv , where 
L is given by (3.15). 

Proposition 3.2. Let q G L°°(0), let (3 be given by (3.16), where (3 G C 4 (0; R+) fulfills 
Assumption 3.1, let if and r\ be as in (3.17), and let L, M\ and M 2 be defined by (3.15)- 
(3.18). Then there exist A > 0, s > 0, and a constant C > depending only on T, 
\\q\\l°° (o) > O an d To? such that the estimate 

s\\e s?? Vu;||l2((_ TiT ) xC ,) + s 3 ||e sr? ^||L2((_r,T)xo) + H^' e SVw \\l 2 ((-t,t)xO) 

J'=l,2 

< C ( S ||e- s V 1/2 (^/3) 1/2 a^||^ (( _ TiT)xro) + \\e-^Lw\\l {{ _ TtT)x0) ) , (3.19) 



holds for all s > sq and all w G L 2 (— T, T; Hg(C?)) satisfying Lw G L 2 ((— T, T) x 0) and 
^6L 2 (-T,T; L 2 (r )). 

3.2 Proof of Theorem 1.1 

In this section we apply the Bugkhgeim-Klibanov method defined in [3] to prove the 
stability inequality (1.6). To this purpose we first introduce the following useful notation 
used throughout the section. For all d > we define := uj x (— d, d). 

Fix L > £, set r = (L + £)/2 and consider a domain O in R n , with C 2 boundary dO, 
obeying 

sit c n r c o c n L . (3.20) 

Next choose x e C^QR^, [0, 1]) such that 





if |x n | 


{o 


if |x n | 



Put io = X v, where v is the C°([-T, T]; Hj(fi) n H 2 (fi)) n C\[-T, T]; L 2 (fi))-solution 
to (2.14), in such a way that the restriction u>i<p of it) to O satisfies 

w\ G C°{[-T, T}- Rl(0) n H 2 (G)) n C 1 ([— T, T]; L 2 (C)) (3.21) 

and is solution to the following system 

-iw' - Aw + giio = xf - in T) x £>, 

u>(0,x) = ix(g 2 - Oi)(x)u (x), x G 0. (3.22) 
w(t,x) = 0, (t,x) G (— T, T) x <9£>, 

where = [A,x] = x + 2x^» n - Here x (resp. x) is a shorthand for the first (resp. 
second) derivative of x- 

As a preamble to the proof of the stability inequality stated in Theorem 1.1 we first 
establish two elementary technical results. 

3.2.1 Two auxiliary results 

Lemma 3.3. For s > 0, let <fi := e~ sr, w, where w is defined by (3. 21) -(3. 22). Then we 
have 

J ■= We'i^wiO, = 21m ( / M^dtdx] . 

Proof. In light of (3.16)-(3.17) we have lim r)(t,x) = +oo for all x G O hence 

lim <b(t, x) = 0. 

tl(-T) 
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Therefore J = ||0(O, OIIl^O) = L-to)xO dt\4>\ 2 dtdx, from where we get that 



J = 2Re ( / d t <p<pdtdx ) . (3.23) 

'(-T,0)xO 



On the other hand, (3.18) and the Green formula yield 
Im ( / (Mxcfy^dtdx) 

\J(-T,0)xO J 

= Re ( / d t (f4>dtdx ) + Im ( / A<jxj)dtdx + s 2 \\ V77<^||l2(_ Ti0)xC , ) 

\J(-T,0)xO J \J(-T,0)xO ' J 



'(-T,0)xO 

I ' , w 

'(-T,0)xO / \J(-T,0)xO 

Re I / d t <jxf)dtdx I + Im ( || V0||l2((_ t ,o)xO) ) = ^ e ( / d t 4>4>dtdx I , 



so the result follows from this and (3.23). □ 
Lemma 3.4. Let w and J be the same as in Lemma 3.3. Then we have 

J < s~ 3/2 I(w), s > 0, 

where 

I(w) :— s\\e sr? Vu>||L2((_ TiT ) xC ,) + s 3 ||e s ^HIl 2 ((-t,t)xO) + H^j e SVw Wl 2 {(-t,t)xO)- 

J"=l,2 

(3.24) 

Proof. In view of Lemma 3.3 and the Cauchy-Schwarz inequality, we have 

J < 2||Mi0|| L 2((_t,t)xO)II0IIl 2 ((-t,t)xO) 

< s~ 3/2 (||M 1 (e- s ^)||^ (( _ TiT)x0) + s 3 ||e^||2 2(( _ TT)x0) ) . 

This and (3.24) yields the desired result. □ 
3.2.2 Completion of the proof 

The next step of the proof involves majorizing I(w) with the aid of the Carleman 
inequality of Proposition 3.2. In view of (3.21)-(3.22) and (3.24) the estimate 



I(w) < C[s\\e 



- s \ l, \d u ^l 2 d v w\\h {{ -T^xT 0) + He^V - Kv)\\l {{ _ TtT)x0) ) 



< C (s\\e~ s \y 2 (d u /3)y 2 dM\h (i -T,T)xr o) 



+ ^2 H e SV ^ j v \\L 2 ((-T,T)x(n r \n e )) + ll e SV xf'\\h 2 ((-T,T)xn r ) I > (3.25) 

j=0,l J 



holds for s > sq. Here and henceforth C denotes a generic positive constant. Bearing 
in mind that rj(0,x) = mfte(-T,T) n(t,x) for all x G Q and that /' = (q 2 — qi)u' 2 , u' 2 , we 
have 

lle-^x/'H^cc-T^xQ.) < C\\e^\ qi - g 2 )|| LW (3.26) 

Here we used the fact that ||w 2 ||L°°((-T,T)xn<!) < 00. On the other hand, the function 
<pd v {3 being bounded in (— T, T) x T as well, we deduce from (3.25)-(3.26) that 



I(w) < C ( sWe-^dMlha-T^ro) + E lle-^^-V^' 

V 3=0,1 



2 

t,|l L 2 ((-T,T)x(O r \Q^)) 



+||e-^( gi - 9 2)||^ (a ,) 

< c(obs + ||e-"(°'% -92)11^)) , s>s Q , (3.27) 

where 

obs := S ||e-^°'^^||^ (( _ TiT)xro) + £ l|e-^ (Ol0 V J "t;||^ (( _ Tir)x(npXfl|)) . (3.28) 

3=0,1 

In light of Lemma 3.4, (3.27)-(3.28) then yields 

J = \\e- s ^\( qi - q 2 )u \\h {o) < Cs^ 2 (obs + ||e-"<°-% - q 2 )\\h m ) , (3.29) 

for all s > s . Further we have J > a 2 \\e~ sri ^ ''\qi — q 2 ) 1 1 l 2 (o £ ) fr° m the assumption (1.3) 
whence 

(a 2 - Cs- 3 / 2 )\\e- s ^\ qi - q 2 )\\ 2 L2m < Cs'^obs, s > s , (3.30) 

by (3.29). Chosing s > Sq so large that a 2 — Cs~ 3 ^ 2 > a 2 /2 and taking into account 
that mi x&ne e- sr > {0 ' x) > 0, we derive from (3.28) and (3.30) that 

\\Qi ~ Qa\\iP(si t ) < C \\d„w\\h((- T ,T)xr ) + Yl II V ^Il 2 ((-t,t)x(ivW) I • ( 3 - 31 ) 
V 3=0,1 / 

Now, recalling from §2.3 that \\V j v ||l2((-t,t)x(qa^)) = 2 II V J 'u||La(( ,T)x(n P \n < )) for 3 = 
0,1, and from the identity w = x v that ||9 J/ w||L2((_T i r)xr ) — 2||<9^u>||L2((o,T)xr ), it 
follows from (3.31) that 

lift - gallon) ^ C \\ d M\h((o,T)xr ) + E II V ^Hl 2 ((o,t)xOM^)) I ■ ( 3 - 32 ) 

V 3=0,1 / 

Here we used the identity \\qi - 92 Hl^) = ||?i - feHi^n) arising from (1.4). 
Set Ti := du x (— r, r) and T 2 := T R (fix \ fi r ) in such a way that 

r = riur 2 , rinr 2 = 0, 
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in virtue of (3.20). Since x( x n) — for every \x n \ > r then the function w = \v is 
uniformy zero in (0, T) x (JIl \ &r), entailing: 

d u w(t,a) = 0, (t,a) G (0,T) x T 2 . 

As a consequence we have 

ll^u>||L 2 ((o,T)xr ) = ||^w|| L 2(( 0T ) xri ) < C||9„u|| L 2(( ,T)xri)- 
Putting this together with (3.32) and the identity v — {u\ — U2)' we obtain (1.6). 
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